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Abstract 

The one-dimensional Coulomb-like potential with a real coupling constant [3, and a 

2 1 

centrifugal- like core of strength G = a 2 — \, viz., V(x) = t^-ic) 2 + \x-ic\ ' ls discussed 
in the framework of PT-symmetry. The PT-invariant exactly solvable model so formed, 
is found to admit a double set of real and discrete energies, numbered by a quasi-parity 
q = ±1. 
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I. Introduction 



Ever since it was conjectured by Besis, Bender and Boettcher, and others [1-3] that Her- 
miticity of the Hamiltonian is not essential for the reality of the spectrum, non-Hermitian 
Hamiltonians have attracted a lot of attention. The main reason for this is that a large 
number of one- dimensional complex potentials, invariant under the simultaneous actions 
of space- and time-reflection operators V and T respectively, have been found to admit 
real and discrete energies. Though this unusual behaviour was at first thought to be 
because of the so-called PT-invariance, this condition is neither necessary nor sufficient 
to ensure the reality of the spectrum. Some authors have pointed out that the necessary, 
but not sufficient, condition for the spectrum to be real and discrete is the 77-pseudo- 
Hermiticity, rfKvj' 1 = 7i + , where 77 is a Hermitian linear automorphism [4-7]. However, 
no general condition has been found for the breakdown of PT-symmetry either. Several 
authors have studied many one-dimensional non-Hermitian models, constrained by their 
PT-invariance [2-12], and shown that such Hamiltonians exhibit 2 types of behaviour - 

(i) In the unbroken PT-symmetry phase, the eigenfunctions of the Hamiltonian 7i are 
also eigenfunctions of PT, and the energy spectrum is real and discrete. 

(ii) When PT-symmetry is spontaneously broken, though the potential retains PT- 
symmetry, the corresponding wavefunctions cease to be the eigenfunctions of the operator 
PT, and the energy eigenvalues are arranged as complex conjugate pairs. 

The interesting fact observed about PT-symmetry is that the non-Hermitian, PT- 
invariant models admit some of the properties of the usual Hermitian ones, viz., super- 
symmetry, potential algebra, quasi-solvability, etc. This has motivated many authors to 
study such non-Hermitian Hamiltonians, especially because of their applications in many 
areas of theoretical physics — nuclear physics, quantum field theories, scattering prob- 
lems, localization-delocalization transitions in superconductors, defraction of atoms by 
standing light waves, as well as the study of solitons on a complex Toda lattice [13]. 

In this note we study the PT-symmetric one-dimensional Coulomb-like potential, with 
a centrifugal-like core, at a real coupling (3 : 



The motivation for investigating such a system arises from the fact that the Harmonic os- 
cillator and the Coulomb potentials are the most widely studied of all quantum mechanical 
systems. Though the PT-symmetric version of the one-dimensional Harmonic oscillator 
was explored by Znojil [14] three years ago, the one-dimensional analogue of the Coulomb 
potential is yet to be investigated. It is worth mentioning here that Znojil has regularized 
the Coulomb potential at an imaginary coupling, and developed it from the Harmonic os- 
cillator by the PT-symmetric Kustaanheimo-Steifel (KS) transformation [15]. However, 
the KS-type mappings also change the dimensions and angular momenta and the energies 
of one system are related to the coupling constants of the other and vice versa. Thus the 
Harmonic oscillator in D dimension is mapped into the (i-dimensional Coulomb problem. 




(1) 
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The approach used in the present work is totally different. The Coulomb potential is 
treated in one dimension, with a real coupling constant. 

The organization of the study is as follows. To make the work self-contained, a brief 
review of PT-symmetry is given in Section II. In Section III, the one-dimensional Coulomb 
problem is developed in the framework of PT-symmetric quantum mechanics. The wave 
functions are normalized in Section IV. Section V is kept for conclusions and discussions. 



II. PT- symmetry 

Before proceeding further, let us briefly recapitulate the idea of PT-symmetry. The 
Hamiltonian H for a particle of mass m, in a complex potential V(x) = Vr(x) + iVj{x) is 
given by 

7i is said to be PT-symmetric when 

VTH = HVT (3) 

Here V is the Parity operator acting as spatial reflection, and T stands for Time Reversal, 
acting as the complex conjugation operator. Their action on the position and momentum 
operators are given by : 

P : x — > —x, p — » —p 
T : x — > x, p — > —p, i — > —i 
Hence, in explicit form, the condition for a potential to be PT-symmetric is 

[V(-x)Y = V(x) (4) 

so that Vr(x) must be an even function of x, whereas Vr(x) has to be odd. The commu- 
tation relation 

[x,p\ = ih (5) 

remains invariant under VT for both real as well as complex x and p. 

For the probabilistic interpretation of the wave function of such a non-Hermitian quan- 
tum mechanical system, the norm and the inner product have to be redefined [16,17]. The 
counterpart of the scalar product of two eigenfunctions ijji(x) and ip2(x) is defined as 

<H^i> = J dx ij>Z(-x) fa(x) (6) 
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and the normalization condition is replaced by 

/oo 
dx ip*(-x) ip(x) = a , a = ±1 (7) 

-oo 

It is evident that the norm is not necessarily positive-definite. So it is referred to as the 
pseudo-norm, a = 1 corresponds to the PT-symmetric phase while a = — 1 corresponds 
to the PT-antisymmetric phase. 



III. One-dimensional Coulomb-like potential 

We start with the radial Coulomb problem in three dimensions. The Schrodinger equation 
for this system (in units h = 2m = 1) is 

_ ii 



Changing its variable r to x = r + ic, c > 0, and replacing the angular momentum I by 
a — \ [14], we obatin the PT-invariant one-dimensional Coulomb-like potential, with a 

„2 i 



2 

centrifugal-like core of strength G = a 2 ; 



a 2 - i 8 
V (x) = 7 r% + (9) 

Here (3 is a real coupling constant and ct > 0. The analytic continuation of the model to 
one dimension makes the boundary condition at the origin redundant. The shift c avoids 
the singularity from the integration path, rendering the Hamiltonian exactly solvable on 
the entire line for any value of a > 0. The Schrodinger equation (8) then takes the form 

ft + [ E - Q2 ~i _ t ) - o (10) 
dx 2 \ (x — ic) 2 \x — ic\ J 

After some straightforward algebra, the solution to eq. (10) is obtained as 

V>(*) = Kn^ lX ~ iCl HX - ic]}-** 1 * L-^^X-icl} (11) 

where N^ n is the normalization coefficient 

= KnK (12) 



7 = 2vK Q J (13) 
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and L n 2qa are the associated Laguerre polynomial [18] 

(14) 

l$(x) = 1 (15) 

I™(x) = m + 1 - x (16) 

From eq. (11), ip*(—x) = ift(x), so that this is a case of unbroken PT-symmetry. As 
expected, the spectrum turns out to be real and discrete. The eigenenergies are obtained 
as a double series, parametrized by a quantum number q = ±1, given by 



(2n - 2qa + 1) 



q may be called the quasi-parity of the system. In the limit c — * 0, the Hamiltonian 
gets back its hermiticity, reverting back to the ordinary radial Coulomb problem, and the 
quasi-parity degenerates to the ordinary parity. 



IV. Normalization of the wave function 

In PT-symmetric quantum mechanics, vj)(x) and ip*(—x) both satisfy the time-independent 
Schrodinger equation 

" ^ + V(x)^ = Ei, (18) 

If E is real and non-degenerate, then 

ij>*(-x) = e^tj}{x) , < < 2tt (19) 

For the phase factor = 0, ip(x) is PT-symmetric, and for = it, ip(x) is VT- 
antisymmetric . 

For the particular potential we are dealing with, = 0, so that the functions ift(x) are 
PT-symmetric. Therefore the normalization condition reads 

/oo 
dx ip*(-x) if)(x) = 1 (20) 
-oo 

/oo 
dx p{-x) 4(x) = \N« n \ 2 r qn (21) 
-oo 



Now 
where 



C = r dx e-^-l { 7 |x-^|}- 29Q+1 {L-^ ^\x - ic\)} 2 (22) 

J — oo k ' 



is evaluated by performing the contour integration / dz f(z) in terms of the complex 
variable z — \x — ic\, with f(z) being given by 

f(z) = ["dze-r { 1Z y w {L-^[ 1Z \} 2 (23) 



Im z 




Re z 



Fig. 1. Contour T in the complex plane 



Thus f(z) is analytic everywhere except at the branch point z — 0. Enclosing the contour 
T = PQRSTU in the complex plane, avoiding the singularity at z = by a semicircle of 
radius e, one obtains j r dz f(z) = from Cauchy's theorem. Detailed calculations show 
that in the limit A — > oo and e — > 0, Pf^, f(z) — > 2 J °° /(x). Also, converges 
for qa < 1. Thus for g = ±1, a should be restricted in the range < a < 1. 
After some algebra, is calculated to be [18] 



r 

qn 



2(-2ga + l) T(-2ga + n+l) 
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71! 



1 - h) F {-qa + 1, -ga + |, - 2ga + 1, ^^j) 



(24) 



(1 + /i)-2 9 a+2 

- ^ h=0 

where F are the generalised Hypergeometric functions [18]. Thus the normalization coef- 
ficient is evaluated from 



\n: 



qn I 



qn 



(25) 



It is easy to observe that for the ground state , 



l JV g0l 



T(-2qa + 2) 



(26) 
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and for the first excited state, 



pa 



\(-2qa + 3) T(-2qa + 2)< (2?) 



V. Conclusions and discussions 

The VT- invariant one- dimensional Coulomb-like potential , _ ic , , together with a cen- 
trifugal core ? x ZiJ\i ; has been solved exactly at a real coupling constant [3. In this respect 
the present analysis is different from the imaginary coupling regularization of ref. [15]. 
The shift c avoids the singularity from the integration path. In the limit c — > 0, the 
Hamiltonian reduces to that of the conventional Coulomb potential, and the quasi-parity 
degenerates to the ordinary parity. 

It is observed that, analogous to the PT-invariant harmonic oscillator [14], the eigen- 
functions as well as the eigenenergies are numbered by a quasi-parity q — ±1, and the 
spectrum is found to be real and bounded. Thus the PT-symmetric version gives rise to 
a second set of bound states. Moreover, these energies are found to depend on the value 
of a. a = is not allowed as the quasi-parity no longer exists in this case. Normalization 
of the wave functions applying the modified normalization for 'PT-symmetry, restricts 
the value qa < 1, which for q = ±1, implies < a < 1. While q = — 1 reproduces the 
ordinary three-dimensional Coulomb energies, viz., 

E" = j , n = 0, 1,2,- •• (28) 

(2n + 2a + l) 2 ' 1 ; 

q = +1 sector further enriches the spectrum. For this particular case, a — | is not allowed 
for n = 0, as then the energy disappears from the system. Therefore, 

(2n-2a+iy 



K + = -~ „ , —2 , n=l,2,3,... (29) 



In a way this may be compared with supersymmetric quantum mechanics, where the 
partner potentials are isospectral, with the possible exception of the ground state. 

To conclude, the analysis in the present work is different from the imaginary coupling 
regularization in ref. [15], where one encounters unavoidable level crossings at all positive 
integral and half-integral values of a. Moreover, for q = +1, there are a plethora of flown- 
away energies at a — n-\-\. However, it is shown in this study that a is restricted to lie 
in the range < a < 1, and hence such phenomena are not observed here. The single 
flown-away energy at a = \ for q = +1 can be avoided if one restricts n to n — 1, 2, 3, • • •. 
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